Modified version 
February 1, 2008 



Fermionic determinant as an overlap between 

bosonic vacua. 

CD. Fosco a * 
and 

F.A. Schaposnik 6 ^ 

a Centro Atomico Bariloche, 84 00 Bariloche, Argentina 
b Departamento de Fisica, Universidad Nacional de La Plata 
C.C. 67, 1900 La Plata, Argentina 



Abstract 

We find a representation for the determinant of a Dirac operator 
in an even number D = 2n of Euclidean dimensions as an overlap be- 
tween two different vacua, each one corresponding to a bosonic theory 
with a quadratic action in 2n + 1 dimensions, with identical kinetic 
terms, but differing in their mass terms. This resembles the overlap 
representation of a fermionic determinant (although bosonic fields are 
used here). This representation may find applications to lattice field 
theory, as an alternative to other bosonized representations for Dirac 
determinants already proposed. 
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Based on an earlier idea of Kaplan PL the overlap formalism [0, |3[ has 
been proposed as a way to define fermionic chiral determinants. Its lattice 
implementation seems to overcome the kinematical constraints imposed by 
the Nielsen-Ninomiya theorem [|J . It may thus provide a suitable framework 
to study interesting non-perturbative phenomena in models containing chiral 
fermions. 

In the overlap approach, the determinant of a chiral Dirac operator in 
D = 2d dimensions is defined as the overlap, i.e., scalar product, between the 
Dirac vacuum states of two auxiliary Hamiltonians acting on Dirac fermions 
in 2d + 1 dimensions. This method has been recently extended to the case 
of Dirac determinants in odd dimensions ||. 

The overlap between two vacua can be implemented in at least two dif- 
ferent ways: an 'operatorial' version, based on building up the Dirac vacua 
by occupying all the respective negative energy states, and then computing 
the scalar product; and also in a path integral approach, which introduces an 
extra dimension into the game. The ovelap is obtained by calculating a path 
integral amplitude for a system whose Lagrangian has a mass that depends 
on the coordinate labeling points in the extra dimension (in the spirit of the 
domain- wall picture). 

In this letter we shall provide an alternative definition for a fermionic de- 
terminant as an overlap between two bosonic vacua. The issue of bosonizing 
a Dirac operator, namely, writing a fermionic determinant in D dimensions 
as a functional integral over purely bosonic fieldsf] has recently received aten- 
tion |7], |j. The approaches 0, || share the property of involving an infinite 
number of bosonic fields, although for different reasons than in the overlap 
approach. The infinite number of fields manifests itself as an extra dimen- 
sion in ||, and as a (discrete) infinite number of fields in the approach of [/]]. 
Besides, they deal not with the chiral case, but rather with systems involving 
both chiralities, i.e., Dirac fermions, the same case we shall consider in the 
present work. This is, in practice, the interesting case for lattice QCD. In 
treatments that avoid the quenched approximation, a bosonic representation 
for the quark determinant (which is real) may be very useful indeed. We 
will, in this letter, also deal with this case. 



1 This is all what is meant by 'bosonization' here. Note that this meaning is quite 
different to the one used, for example, in ||. 
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In our construction, it is convenient to write the vacuum state in the func- 
tional Schrodinger representation [10| (we shall later on discuss an operatorial 



representation). In the Schrodinger representation, the " wavefunction" that 
describes the vacuum state for a real scalar field tp with a quadratic action 
is, in the "coordinate" basis, a vacuum functional (<f\^o) = ^/q (if), of the 
following kind 

*n(v) = det^(-) exp{-~(pQ<p) (1) 

7T I 

where we are using a shorthand notation (similar to the one of [fL2|), such 
that 

<pSl<p = J d D xd D y<f(x)n(x,y)tp(y) (2) 

and Q(x, y) is a real, symmetric, and definite positive kernel. 

In the coordinate basis, all the wavefunctions depend on the field configu- 
ration <p(x), where x labels the D 'spatial' coordinates in a D + l dimensional 
spacetime. Thus the scalar product between two states $1, $2 is defined by 
a functional integral 

($i|$ 2 > = / VtpQXip)^) (3) 

where the integration measure in (|3|) is formally defined as Dtp = Yl x dip x 
with x in the D dimensional space. 

The functional determinant factor in (|I|) is introduced in order to nor- 
malize this wavefunction 

(*n|*n) = 1 • (4) 

This vacuum state can be thought of as the ground state for a (second 
quantized) Hamiltonian if) 



d u x7r 2 (x) + / d u x / d u yif(x)Q 2 (x,y)if(y) 



(5) 



with 7r(x) = 8/i8f{x) the momentum operator, and f acting multiplicatively. 
Following notation (||), we write Q 2 (x,y) = J d D zQ(x, z)fl(z,y). As an 
example, one has for a free real scalar field Q 2 (x,y) = (— V^, + m 2 )5(x — y). 

Let us now consider the overlap between two vacuum states, ^^(y?), 
tyn 2 (tp), corresponding to two Hamiltonians, differing just in their (com- 
muting) kernels Q\, ^2- Performing a functional Gaussian integration, and 
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rearranging terms, the overlap (^q 1 \^n 2 ) yields 



^nj*^ = det3(^)det-5(^— 



det 4 



1,^1 ^ 2 

4<^ + s; +2 » 



(6) 



Or, by defining O = 

(* ni |*n 2 ) = det* [/(£>)] (7) 

where 

/( *> - jrhr • (8) 

and we have ignored (and so will do in what follows) irrelevant constant 
factors. Determinants are supposed to be adequately regularized 

so that all usual properties (like det(AB) = det A. det B) ) hold. The 
symmetry f(x) = /(-) is a consequence of the symmetry of the overlap (|^) 
under 1^2. 

The generalization to the case of a complex scalar field, endowed with 
a complex hermitean definite positive covariance Q, is straightforward. One 
obtains, instead of ([7p, the result 

(* ni |*Q 2 ) = det*[/(0)] . (9) 

This of course generalizes to the case of Nf complex flavours. The results 
implies just to a change of the exponent of the determinant in the form (0) 

N f 

todet-[/((9)]. 

As usual in similar approaches, we shall assume that the operators Q\, 
2 have been regularized (on the lattice, for example), what makes them 
bounded. We shall also assume there is a non-vanishing gap in their spectra. 
Then their ratio O is also regularized. Moreover, we will let it depend on a 
mass parameter, say M, such that 

M -»• oo \\0(M)\\ << 1 (10) 

then 

(*nj*n 2 ) ~ det5[-l] when M -> oo . (11) 

1'2 
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To deal with the specific example of Dirac operators, we will consider the 
case of fli and Sl 2 such that 



nj = -p 2 + m 2 , n 2 2 = -I/> 2 + M 2 (12) 

where p = = + A^, and A^ is an (Abelian or non-Abelian) 

gauge connection. To apply the result flllD to the present case, we need to 
assume that is regularized (in the lattice, say). Then, keeping the regulator 
cutoff finite, and letting M >> m, M >> \ \ J/>\\, we have, for the case of Nf 
flavours, 

/t it x i N f r -$ 2 + m 2 1 , N f ,-B 2 + m 2 , . . 

(*n, *na> ~ det * — S ~ det* ^ J . 13 

Note that the regulating cutoff must be kept finite when taking this large M 
limit, and indeed, M should be in fact bigger than the cutoff. Of course, if 
we normalize by dividing ( |T3"D by the same object evaluated for zero external 
field (A = 0), we obtain 



(*n,l*(i 3 )U.o '-^ + m 2 ' 

Of course, the number Nf of bosonic flavours in the overlap must be adjusted, 
depending on the actual number of fermionic flavours Np considered. It is 
evident that we can only consider an even Nf, and moreover Nf = 2Np. 
Then,we shall write our main result as 

det [ ^ 2 ] ~ /xT , ,. Tf ri (15) 

We have then been able to write the Dirac operator determinant for a D- 
dimensional theory of Dirac fermions with Np flavour in terms of the ratio of 
two overlaps between vacuum states in a D + 1 dimensional bosonic theory. 
Note that since we are working with a Dirac operator in D = 2n dimensions, 
one has the identity det(— lp 2 + m 2 ) = det(^ + m) 2 . Then, our overlap 
formula gives in fact the value for \det(ft + m)\. We see that our approach, 
as Slavnov's and Luscher's, only works in even dimensions. The reason is 
that only in even dimensions is the Dirac determinant real. The extension of 
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our method to odd dimensions would imply the necessity of defining vacuum 
functionals for non-Hermitian Hamiltonians. 

There is a striking similarity between this case, for a finite M (see equa- 
tion (^)), and what the standard overlap yields for the modulus of a chiral 
determinant: \(A + \A—)\ = det 5 ( C ^+ M2 1 , where C is the chiral Dirac 
operator [fT3| . 

For the concrete example of the Dirac operator in D = 2 dimensions, 
we may even interpret the ratio O as a Pauli-Villars regularization for the 

operator ../Oi = \J — 1J) 2 + m 2 , since just one regulator field suffices (the 
contribution of this regulator is of course f2 2 )- We still assume a lattice 
regularization. Then the overlap for this case is doubly regulated, when one 
takes the large-M limit (in particular, M >> 1/a, where a is the lattice 
spacing), one is removing the Pauli-Villars regulator, and at the same time 
approaching the determinant. The latter is still regulated since a is kept 
finite. 

We will also give an operatorial construction of the overlap (||) between 
the two vacuum states, which avoids the use of the Schrodinger representa- 
tion. We need to introduce suitable creation and annihilation operators. As 
the Hamiltonian is quadratic, we shall first consider the case of one simple 
harmonic oscillator, and then generalize to the case of interest, which contains 
an infinite collection of decoupled harmonic modes. For a harmonic oscilla- 
tor, the analogous of @ would be to evaluate the scalar product between 
to vacua, each one corresponding to a Hamiltonian with a given frequency 
u. Namely, we need to evaluate (using operatorial methods) the object 
(tyui\^u> 2 ) where \^f ui ) and l^^) are the ground states of Hi = |(p 2 + oj\q 2 ) 
and H 2 = \{p 2 + uj 2 q 2 ), respectively. It is possible to show that, if we de- 
fine the corresponding two sets of creation and annihilation operators Oj, aj, 
% = 1,2, they will be related by the Bogoliubov transformation. This trans- 
formation is exactly of the kind that appears when defining " squeezed states" 

di = cosh6*a 2 + sinh^aij = Uga 2 Ug 

aj = sinh6 l a2 + cosh^a^ = Uga 2 Ug (16) 

where 



cosh 6 = 
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^ - iWS-^J (17) 

and 

U e = exp[^(a 2 - (4) 2 )] . (18) 

^From this it follows that we can write one vacuum in terms of the other as 
^ui) — Uel^/^). It is evident that I^J, when represented in the Hilbert 
space built with a\ acting on |^ 2 ), will be a linear combination of states 
containing an even number of excitations. After some algebra, this can be 



put more explicitly as follows [13 



--In cosh 9 --(a^) 2 tanh# 
= c 2 e 2 K J |^ 2 ). (19) 

^From (FTO]), we obtain for the overlap 

--In cosh 9 oj-t x 

(*J*W = e 2 = [ 2/ (-l)]4 (20) 

with / as defined in (|j). This is clearly the equivalent of (|7J) for the case 
of a + 1 field theory. But the generalization to the D + 1 case is trivial 
because the Hamiltonians are quadratic, and the system is brought to an 
infinite collection of uncoupled harmonic oscillators by using operators that 
create or destroy particles occupying states that are eigenmodes of the kernels 
Q 2 and Q 2 . (they are Hermitian, and differ by an operator proportional to 
the identity). These, of course, are the usual plane waves when Q 2 (x,y) = 
(-V 2 , + m 2 )5(x - y). 

Defining thus the eigenmodes g k {x) by 

d D ytt 2 1 (x,y)g k (y) = \ 2 k g k (x) , d D ytt 2 2 (x,y)g k (y) = \' 2 k g k (x) (21) 



we obviously have \\ — A' 2 , = M 2 — m 2 . Then we can apply (pOf ) to each 
mode, which now has frequency Af. (or A' 2 .). This yields 

(^nj^) = IMp)] 4 = deti[/(0)] (22) 
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in agreement with (0). 

We conclude mentioning that our approach shares with the overlap def- 
inition and with the proposals of || 0] the property of involving an infinite 
number of fields, here manifested in the use of a system living in D + 1 di- 
mensions. Our method uses bosons rather than fermions, what means that 
the D + l dimensional theory would be rather unphysical, as in 0, f|, in the 
sense that the action attributed to the bosons would be quite exotic. It may 
however, be a useful technique for lattice simulations, even in the framework 
of the Schrodinger functional approach, which can be implemented on the 
lattice [[□]] . This because our representation would allow to include the effect 
of fermion loops without having to simulate fermions, a well known stumbling 
block for numerical simulations. Of course, one may also attempt to solve 
for the vacuum functional analyticaly, for example by variational techniques, 
and then calculate the overlap. 
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